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Abstract 

We argue that the integrable modified nonlinear Schrodinger equation with the 
nonlinearity dispersion term is the true starting point to analytically describe 
subpicosecond pulse dynamics in monomode fibers. Contrary to the known 
assertions, solitons of this equation are free of self-steepening and the breather 
formation is possible. 

1. Introduction 

Soliton-based optical communication systems serve as an exciting example 
of the application of a purely mathematical concept (soliton) to modern tech- 
nology. The nonlinear Schrodinger equation (NSE) 

iu t + (l/2)u xx + \u\ 2 u = (1) 

is the adequate model to describe picosecond soliton evolution in monomode 
fibers 0| . Here u(x, t) is the envelope of the pulse electric field and coordinates 
t and x measure distance along the fiber and time in a frame comoving with the 
pulse group velocity, respectively. The applicability of NSE depends crucially 
on the assumption that the spatial width of the envelope is much larger that the 
carrier wavelength. Besides, the success of this model is substantially related to 
integrability of NSE || and hence to the controllability of soliton parameters § . 
Various more subtle effects accompanying the picosecond soliton propagation are 
usually treated as a perturbation of the integrable model. 

On the other hand, dynamics of subpicosecond optical pulses (< 100 fs) is 
not well governed by NSE because the above mentioned assumption is not sat- 
isfied. The spectral width of subpicosecond pulses becomes comparable with 
the carrier frequency, and three main additional effects - nonlinearity disper- 
sion, intrapulse Raman stimulated scattering and linear higher-order dispersion 
- should be taken into account [Q: 

iu t + (l/2)u xx + \u\ 2 u = iai (\u\ 2 u) x + a 2 (\u\ 2 ) x u + ia 3 u xxx . (2) 

The terms in rhs of (2) account for the above additional effects. In general, extra 
terms violate integrability of the equation. Hence, a question can be posed: does 
there exist an equation that will be integrable as NSE and at the same time 
would be more relevant in the subpicosecond range? The answer is positive 
because the modified NSE (MNSE) 

iu t + {l/2)u xx + \u\ 2 u + ia (\u\ 2 u) = 0, a e Re (3) 
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with the a-dependent nonlinearity dispersion term is still integrable though 
the associated spectral problem is more involved than the Zakharov-Shabat 
one. Namely, the initial-value problem for MNSE (3) can be solved within 
the framework of the Wadati-Konno-Ichikawa (WKI) spectral problem 0. A 
careful study of the WKI spectral problem (or the quadratic bundle) for MNSE 
and related equations was undertaken by Gerdjikov and Ivanov Q. Explicit 
soliton solutions to MNSE obtained in and ^ turned out too complicated 
for practical use. That is the reason that MNSE is usually treated as NSE with 
a-dependent perturbation term, especially as actual values of a are normally 
small. 

Among the other things, such treatment gave rise to some misunderstand- 
ing. First, it was shown in j^] that the initially symmetric hyperbolic-secant 
pulse evolving in accordance with MNSE (3) develops an asymmetric self-phase 
modulation and a self-steepening. There is a wide-spread opinion that the self- 
steepening is an inherent property of the subpicosecond pulse dynamics [|| that 
should be minimized for proper operation of an information system. Second, 
it is well known that the initial pulse 2sechx evolving according to NSE 
produces the NSE breather (the bound state of two solitons). On the other 
hand, the same initial pulse decays into separate solitons when evolving accord- 
ing to MNSE. Hence, it was inferred that MNSE does not admit breathers (or 
higher-order solitons) | Q . 

Our aim here is to show that the situation with subpicosecond soliton dy- 
namics is rather different. We argue that the integrable MNSE is the true 
starting point to analytically describe this dynamics. It is remarkable that nu- 
merical simulations of the MNSE-based soliton propagation revealed various 
regimes which cannot be explained by treating the a-dependent term in (3) 
as a perturbation of NSE [Q. We derive the MNSE soliton solution that is 
non-perturbative w.r.t. a and demonstrate that the MNSE soliton propagates 
without any self-steepening. Besides, we explicitly obtain breather solution to 
MNSE. Numerical simulations confirm the stability of the MNSE breather. 



2. MNSE soliton 



We will employ the Riemann-Hilbert (RH) problem Jl3| for solving nonlinear 
equations. Let us start with the Lax pair for MNSE (3): 

9 x =A(k)[a 3 ,9]+2ikQ9, (4) 

/4i i \ 

9 f = ft(fc) [(7 3 , *] + -k 3 Q + 2tk 2 Q 2 a 3 - -kQ + kQ x a 3 - 2iakQ 3 9, (5) 
\ a a J 

Here the Hermitian matrix Q = f ^ J stands for the potential of the spectral 
problem (4), k is a spectral parameter. The standard procedure is: 
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a) building the Jost solutions of the linear spectral problem (4); 

b) building the solutions <j>± which are analytical in complementary regions of 
the fc-plane; 

c) formulation of the RH problem for Q± with the standard normalization 

$± / for |fc| -> oo, (6) 

where / is the identity matrix. 

It is, however, easy to see by substituting the asymptotic expansion w.r.t. 
to k^ 1 of \f to the spectral problem (4) that this problem does not agree with 
the standard normalization. On the other hand, an associated equation with 
the fifth-order nonlinearity, 

iv t + (l/2)v xx - iav 2 v x + \v\ 2 v + a 2 \v\ 4 v = (7) 

has the Lax pair as well with the WKI spectral problem 

= A(k)[a 3 ,^ A) ] + (2ikQ A + iaQ 2 A as)^ A \ (8) 

- (° S) 

that agrees with the standard normalization, and with the same dispersion 
relation Q,(k) for the temporal Lax equation. Moreover, equations (3) and (6) 
are gauge equivalent and solutions of MNSE (3) follow from those of (6) by 
means of a simple algebraic relation 

Q = 9~ 1 Qa9, g(x,t) = ^ A \k = 0,x,t). (9) 

The associated equation (6) does not have such an obvious physical interpreta- 
tion as the MNSE but it has an extremely simple soliton solution. Hence, we 
will not solve MNSE directly. Instead we will integrate the associated equation 
(6) and then will obtain solutions of MNSE by the algebraic relation (8). 

We begin with the spectral problem (7) for the associated equation. At first 
we define the Jost solutions J±, J± — > I at x — > ±oo, which are interrelated with 
the scattering matrix S, J_E = J + ES. Here E = exp(A(fe)a;). Dividing the 
Jost solutions into columns, J± — (J± \ J±^), it can be shown by the standard 

analysis of integral equations that the columns and are analytical 
in the first and third quadrants of the fc-planc. Hence, the matrix function 
$+ = (J ( +\ J (2) ) is analytical as a whole in the same quadrants. The matrix 
$ + can be expressed in terms of the Jost solution J + and some entries of the 
scattering matrix: 

Because the potential Qa is Hcrmitian, we have the involutions [J±(fe)]* = 
[J±(fc)] , S(k) Jf = S(k)~ 1 . They allow us to introduce the matrix function 
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<&Z 1 (k) = < E>+(fc) t = ES+E^ 1 JZ 1 that is analytical in the second and 
forth quadrants. Thereby, we can pose the RH problem with the standard 
normalization, 

$Z 1 ®+=EG(k)E- 1 , $±->J at fc^oo, (10) 

where 

G(k) = S f + S + = ( _^ S J 2 ) , k E{k = t-i V , Hv = 0} , 

as a problem of analytical factorization of the non-degenerate matrix G(k) de- 
fined on both the real and imaginary axes of the fc-plane. 

In general, the function <J> + has zeros in some points kj lying in the first and 
third quadrants, det <£+(&.,■) = 0. Hence, in these points there exist eigenvectors 
\rij) with zero eigenvalue. It is important that zeros kj appear by pairs (kj, —kj). 
It is a feature of the WKI spectral problem. Hence, the single soliton of the 
associated equation is determined by two zeros k\ and —k\. The zeros kj, 
eigenvectors \n,j) and the matrix G(k) comprise the RH data. Because we deal 
with the solitons only, G(k) being related with the continuous spectrum of the 
spectral problem, is taken to be the identity matrix. 

If <!>+ is a solution of the RH problem (9), it can be expanded in the asymp- 
totic series $ + = I + /& + <fr+ /fc 2 + ■ ■ •. Substituting this expansion into 
the spectral problem (7), we reconstruct the potential Qa'- 

Qa = (l/a)[<7 3 , - (S/aVs^. (11) 

Now we derive a soliton of the associated equation (6). Let us have zeros k\ 
and —hi and two eigenvectors \n±). It can be easily shown that the eigenvectors 
obey the equations 

| a = A(fci)(7 3 |?i + ), \n+) t = f2(fci)<7 3 |n+). 

Hence, we obtain explicit space and time dependencies of the eigenvectors, 

cxp[A(fci)ir + Q(ki)t] exp(a + i<po) 
exp[— A(ki)x — Q(ki)t] 

|n_) = a- 3 |n+), (n±| = \n±) f . 

Here a, <f>o = const. It can be shown by the dressing method jl3) that the matrix 
$+ is represented as (k± = ±fci) 
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Because the eigenvectors are known explicitly, we can evaluate the matrix <I> + 
as $+(fc) = I - D + /{k - h) - D-/(k + fe x ), where 

k\-k\ ( e z (k ie - z + k 1 e z )- 1 e lip {kie- z + k 1 e z )~ 1 



2 \ e-^(fcie 2 + k 1 e- z )- 1 e~ z (k x e z + k^ 2 )- 1 
D- = -a 3 D + a 3 . 

We introduced here new independent variables z and tp: 

z = -(l/w)(x -Vt- x ), <p = Vx~ {1/2){V 2 - w~ 2 )t + <p , 

xq = aw, tpo = const, 
where the soliton velocity V and width w (see below) are defined by 

v = ±(i-2(*; + e;,), — 

Hence, the eigenvalue fci is expressed in terms of velocity and width as 

ki = {l/2){\ - aV -ia/w) 1/2 , Imfci < 0. (13) 

Expanding then $ + in the asymptotic scries, we obtain from Eq. (10) the soliton 
solution of the associated equation: 

(14) 



w kie~ z + k\e z 



It has indeed a very simple form. 

An important aspect of the solution (13) should be noted. Namely, the 
parameter a that enters the soliton width w appears in the denominator. Hence, 
we cannot reproduce the soliton (13) considering the associated equation (6) as 
the a-perturbed NSE. Nevertheless, there exists a procedure |14| to perform the 
limit a — > 0. Namely, representing ki as fci = (1/2) — (a/2)Ai + 0(a 2 ), we 
obtain in this limit from Eq. (13) the NSE soliton with the eigenvalue Ai. 

As regards the MNSE soliton u s , it follows from v s by means of the algebraic 
relation (8). Indeed, u s = (<72/ffi)u s , 

J2 / kl \ D + 22 



Explicitly we have 



ki k\e z + k\e z k\ k%e z + k\e z 



9i = T"i — , , t — 7' 92 



ki k\e z + k\e z k\ k\e z + k\e z ' 

i k 1 e~ z + k 1 e z ilf> 
w (k\e z + kie~ z ) 2 



(15) 
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The MNSE soliton (14) looks much simpler than those derived in |7) and S. 
In the limit a — > 0, the solitons of both the MNSE and the associated equation 
reproduce one and the same NSE soliton. 
Square of module of u s (14) is written as 

K| 2 = \{kie~* + fcie'J-^fcie* + fcie^)" 1 



1 

2^2 



1 - aV + \ (1 - aV) 2 + -J cosh- (a: - Vt) 



We see from this relation that the envelope \u s \ moves holding its shape, i.e., 
without any self-steepening. 



3. MNSE breather 

To derive the MNSE breather, we start from four zeros ±fci and ±k 2 , where 
kj = (l/2)(l — aVj—ia/wj) 1 ' 2 , Ivakj < (cf. Eq. (12)). Because we seek for a 
bound state of two solitons, we put V\=V% = V and, without loss of generality, 
V = 0. We have four eigenvectors with the property \n 2 j) = cr 3 |n 2 j_i), j = 1, 2. 
Namely, 



i exp[—(x/2wi) + i(t/Aw\) + iip n 
"" ~ 1 exp^wO-^t/^+t^ia] 



_ / exp[-(x/2w 2 ) + i(t/4w%)+i(p2i] \ nf ,, 
1 3 > ~ { exp{(x/2w 2 )~i(t/4w 2 2 )+i<p 22 } J [W) 

with a special relation for the constant phases <pn — ipi 2 — <p 2 i+ <p 22 = 7T- Denote 
X 2 j-i = kj and = —kj, j = 1,2. Then the matrix function <&+(&) for the 
breather is written as (cf. Eq. (11)) 

^ n\ t \ n m)(D 1 ) mn (n n \ (n m \n n ) 
$+(k)=I- ^ — , D mn = - 

m,n=l A ™ An m 

We omit cumbersome but evident calculations performed along the lines of Sect. 
2 and give below the explicit expression for the MNSE breather at rest: 



u br = (9'2/9'^Vbr, (17) 
Wl (k ie x/wi + k 1 e' x/wi ) e u / 2w * 



wi - w 2 ! 

v br = ■ D 

Wi + w 2 



D = 



w 2 (k 2 e x l W2 + k 2 e- x,W2 ^j e it/2w * 
Wl w 2 (fcie^/™ 1 + k ie - x / Wl ^j (k 2 e x/w2 + k 2 e- x/v 
w\ (k 1 e x/w +- it/w+w - -k 2 e- x/w++lt/w+w -^j 
( t k 2 e x / w++it / w+w - — k ie - x / w +- it / w + w -^j ) 
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9i 



*+ii(fc = 0) = 1 



2D 



wi - w 2 ( fci / 

1 w iT W ^T~ I e 

wi +w 2 \ k 2 ki 



k 2 kx 

w+ ( ^_ e -2it/w+w- _j_ ^2 e 2it/ u)+ «,_ 

Vfci k 2 



$+22 (fc = 0) = 1 



2D 



'U'l 



+ „■ I !i e 2 «/»+^- 
fcl 



&2 -2it/t 



2x I W-i 



Wl - W 2 ( fcl 

1 w iT w 2i— I e 

Wl + W 2 V fc 2 «1 



2 \ „2a ; / 1U+ 



Here iwjjj 1 = (1/2) (wf 1 iw^ 1 ). It is seen that the MNSE breather oscillates 
with the period T = irw+W- and reproduces in the limit a — > the well known 
NSE breather pOj . Fig. 1 shows the square of module of the breather solution 
(16) for wi — 1/3, w 2 = 1 and a — 0.1. We see that there is no any decay of 
the MNSE breather. 



Conclusion 

We consider MNSE as a natural integrable generalization of NSE to the 
range of subpicosecond optical pulses. It is shown in this paper that MNSE 
possesses the basic ingredients (solitons and breathers) of integrable nonlinear 
equations. To justify the applicability of these results to the description of ac- 
tual subpicosecond pulses, we should account for at least the intrapulse Raman 
scattering that breaks integrability of the equation. A possibility to reduce an 
adverse action of this effect is discussed in |D| on the basis of the perturbation 
theory for the MNSE soliton juj. A novel way to suppress the Gordon-Haus 
effect for the MNSE soliton was revealed in |l6| . Recently quasiradiation solu- 
tion of a compound model including MNSE was obtained by Zabolotskii fl7j| . 
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